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هدف اصلی این پژوهش آشکارسازق رفتار دوره‌ای خشکسالی‌های ماهانه ایران است. برای رسیدن به این 
هدف ابتدا بارش روزانه ۶۱ ایستگاه همدید کشور برای یک دورة آماری ۳۱ ساله (۱۹۸۳-۲۰۱۳) از 
شارشان هراس کین اف کوهه ی ای کال مفن رای هر ای رسای تاهانه 
محاسبه و رفتار دوره‌ای سری‌های زمانی شاخحص خشکسالی مژثر ([(۳1) ایستگاه‌ها با استفاده از تحلیل 
موجک مورد بررسی قرار گرفت. نتایج این تحقیق نشان داد که دوره‌های تناوبی موجود در سری‌های 
زمانی شاخص خشکسالی موثر (۳101) ایستگاه‌های مورد مطالعه در ایران بسیار متنوع بوده و از بازه‌های 
تناوبی ۲۸- ۲۶ ماهه تا ۲۱۰-۲۲۰ ماهه را شامل می‌شود. علاوه بر این. در تعدادی از ایستگاه‌های مورد 
مطالعه» بازه‌های تناوبی غالب به چشم می‌خورد. از جمله بازه‌های تناوبی غالب می‌توان به دو بازه ۱۲۰- 
۸ ماهه و ۸۰- ۰ ماهه اشاره کرد. ولی در مجموع درکل ایستگاه‌ها اين بازه‌های تناوبی غالب مشاهده 
نمی‌شود. همچنین در اکثر ایستگاه‌ها بازه تناوبی در طول سری زمانی نغییر می‌کند. در تمام نواحی به‌جز 
ناحیه نیمه شمالی ایران مرکزی» پهنای نوار تناوبی در حال افزايش است. به عبارت دیگر نوار تناوبی به 
سمت بازشدگی است. در اوایل سری‌های زمانی. نوارهای تناوبی دارای پهنای کمتری بوده و دوره‌های 
کوتاه‌تری را شامل می‌شوند؛ اما با نزدیک شدن به اواخر سری‌های زمانی پهنای نوارهای تناوبی بزرگتر 
شده و دوره‌های بلندتری را شامل می‌شوند. بنابراین در کل ایران می‌توان نتیجه گرفت که بازه‌های تناوبی 
با شیب خیلی ملایمی به سمت دوره‌های بلندمدت متمایل می‌شود. 

کلیدواژه‌ها: خشکسالی, تبدیل موجک. بسامد بازهُ تناوب. 
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۱- مقدمه 


کشور پهناور ایران با بارشی در حدود یک سوم بارش جهانی در یکی از مناطق خشک و نیمه خشک جهان قرار 
گرفته است و خشکی جزهء فطرت و صفت اصلی اقلیم آن به شمار می‌آید. در سال‌های اخیر بنا به عللی که غالبا 
مربوط به تغییرات جهانی اقلیم هستند. ناهنجاری‌های بارش در نواحی مختلف کشور فزونی يافته است و نوسان‌های 
شدید مکانی و زمانی خشکسالی. خسارات هنگفتی را بر اقتصاد کشور به ویژه در طی دهه‌های اخیر تحمیل کرده 
است. کاهش عملکرد مراتع. کاهش تولید محصولات زراعی به ویژه دیم کاهش منابع آب کشاورزی و شرب 
کاهش منابع آب سطحی و زیرزمینی. طغیان آفات و بیماری‌های گیاهی و دامی, افزایش مهاجرت و در نهایت اثرات 
سوء زیست محیطی, اقتصادی و اجتماعی از جمله آثار منفی خحشکسالی‌ها هستند که توسعه پایدار کشور را مورد 
تهدید قرار داده است. 

تغییرپذیری پدیده‌های اقلیمی در بازه‌های زمانی سه‌گانه و به شکل روند. نوسان و افت و خیز خود را نشان می- 
فقتل زرونله رقتار بل یسرم مان و افت دی یر ررفازهای وکا از تخراز مت رانیان ی ند ابا معشته شی او 
کننده نوسان در الگوی تکراری آن است. در این حالت سیستم اقلیمی پس از مدت زمانی پیکربندی "شبیه " به 
الگوی قبلی را با کمی تفاوت به خود می‌گیرد. بسیاری از نوسان‌های اقلیمی در یک دوره ثابت. مثلاً یک سال یا کمتر 
به صورت دوره‌های ماهانه یا فصلی مشخص می‌شوند. اين نوع چرخه‌ها را چرخه‌ها با دوره‌های آشکار می‌نامند. 
بسیاری از چرخه‌های دیگر اقلیمی به صورت آنی به منحصه ظهور نرسیده بلکه به تدریح در طول زمان بروز 
کرده‌اند. این نوع چرخه‌ها را چرخه‌های نهان می‌نامند. 

تاکنون روش‌های آماری و غیر آماری متعددی توسط محققین مختلف جهت شناخت رفتار دوره‌ای متغیرهای 
اقلیمی در قالب مقیاس‌های زمانی مختلف استفاده شده‌اند. یکی از این روش‌ها که اخیراً بسیار مورد توجه محققین 
علوم جوی و منابع آب قرار گرفته است. تحلیل موجک است. استفاده از تحلیل موجک در بررسی‌های اقلیمی و آب 
و هوایی سابقه کوتاهی دارد. اولین کاربرد این روش در هواشناسی را 1060 ک معمصنع (1993) عملی ساختند. 
کاربرد این روش در تحلیل سری زمانی نوسانات جنوبی ال نینو / نوسان جنوبی و محاسبه سطوح اطمینان در تحلیل 
موجک توسط 00۳۳0 ک 10۳6066 (1998) بطور کامل تشریح و مورد استفاده قرار گرفته است. از این تحلیل 
برای توصیف جریان آب رودخانه (1998 ,.21 66 800 رابطه بین شاخحص نوسان جنوبی و مونسون تابستانه 
هندوستان (2000 ,16016810 برای یافتن رابطه بین نوسان اطلس شمالی و تعییرات سطح آب دریا ( ,21 64 ۷۵0 
4 و برای بیان ویژگی‌های سالانه و دهه‌ای بارش (6نط05؟ :2004 ,کموننه ۷ ک دطلقمعط :2003 ,.ل2 اه عصمع‌منک< 
9 ,00عزعمه۹ ک 12/200) استفاده شده است. همچنین ۱۵1620 (1999) برای بررسی تغییرپذیری زمانی بارش 


و رواناب در جنوب ولز و اثر شاحص نوسان جنوبی بر روی آنها از تبدیل موجک پیوسته استفاده کرد. مطالعه این 


محقق نشان داد که یک رابطه قوی آماری بین شاخحص نوسان جنوبی و بارش با یک بسامد غالب ۲۷ ماه وجود دارد. 
صض و همکاران (2000) رابطه بین رواناب و بارش را در چشمه‌های کارستی کشور فرانسه با استفاده از تحلیل 
موجک و تحلیل چند طیفی مورد بررسی قرار دادند. 13011 > 20 (2002) رابطه بین تغییر پذیری تراز سطح 
آب دریا با برخحی متغیرهای هواشناسی را با استفاده از تبدیل موجک گسسته مورد مطالعه قرار دادند و در 
نهایت 16۷0 م2ط) و همکاران (2010) تغییر پذیری پیوند از دور و قابلیت پیش‌بینی بارش را در تایوان با کمک این 
روش بررسی کردند. 

در ایران نیز از تحلیل موجک در شناسایی رفتار دوره‌ای متغیرهای اقلیمی استفاده شده است. تقوی و همکاران 
(۱۳۹۰) با استفاده از این روش تغییرپذیری بارش ماهانه چهار ایستگاه همدید کرمانشاه سنندج. همدان و حرم آباد 
را در غرب ایران مورد بررسی قرار دادند. نتایج کار این محفقین نشان داد که در همه ایستگاه‌ها یک دوره بازگشت 
سالانه با سطح اطمینان بیش از ٩۰‏ درصد وجود دارد. جهانبخش و همکاران (۱۳۹۰) شاخحص دو قطبی دمایی 
اقیانوس هند (10012) و تغییرپذیری بارش‌های فصلی جنوب ايران را با استفاده از تحلیل موجک مورد مطالعه قرار 
دادند. جهانبخش و عدالت دوست (۱۳۸۷) با استفاده از روش آنالیز موجک» نقش لکه‌های خورشیدی بر تغییرات 
بارندگی‌های سالانه ایران را مورد توجه خود قرار دادند. خسروی و میردیلمی (۱۳۹۲) نیز اثر لکه‌های خورشیدی بر 
بارش را در یک مقیاس کوچکتر و به وسعت استان گلستان مورد مطالعه قرار دادند و در نهایت به همان نتایجی 
رسیدند که جهانبخش و عدالت دوست (۱۳۸۷) برای کل ايران رسیده بودند. از تحلیل موجک در تحلیل همبستگی 
بسامدی رطوبت نسبی بابلسر تهران و شاهرود (عراقی و همکاران ۱۳۹۳) و پیش‌بینی خشکسالی‌های ایستگاه 
هواشناسی تبریز (حسن زاده و همکاران» )۱۳٩۱‏ نیز استفاده شده است. 

لذا با بررسی و مطالعه‌ای که بر روی منابع موجود در ادبیات اقلیم‌شناسی و هواشناسی جهان و ایران انجام شد و 
به برخی از آنها نیز در بالا اشاره شده است مشاهده گردید که در بیشتر مطالعات داخلی در استفاده از تحلیل موجک 
همواره کاستی‌هایی وجود داشته است که کمتر به آنها توجه شده است. یکی از این کاستی‌ها مشخص نبودن روش 
تحلیل موجک استفاده شده در این مطالعات بوده است که تنها با اتکا به برخی کلی گویی‌های از ببان آن صرف نظر 
شده است. مورد دوم عدم توجه به این نکته بوده است که تحلیل موجک دوره‌های تناوبی را برآورد نمی‌کند بلکه 
تحلیل موجک بازه‌های تناوبی را برآورد می‌کند. در حالی که در تمامی مطالعات داخلی از این روش برای شناسایی 
دوره‌های تناوبی استفاده شده است. در نهایت سومین نکته اساسی و مهم در تحلیل موجک آن است که وقتی یک 
بازه تناوبی مشخص شد ممکن است این بازه تناوبی در طول دوره وجود نداشته باشد و صرفاً به دوره زمانی خاصی 
مربوط باشد؛ در حالی‌که در مطالعات انجام شده قبلی به این موضوع یعنی پایداری زمانی بازه‌های تناوبی توجه نشده 


و دوره تناوبی برای طول دوره سری زمانی در نظر گرفته شده است. 
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بنابراین با توجه به اينکه تاکنون مطالعه جامع و مستقلی که به طور مشخص هدف آن بررسی و شناسایی رفتار 
دوره‌ای خشکسالی های ایران باشد انجام نگرفته است؛ در این مطالعه قصد بر این است که با در نظر گرفتن 
کاستی‌های قبلی در مطالعات گذشته داخلی اقدام به شناسایی بازه‌های تناوبی خشکسالی های ایران گردد و پایداری 


زمانی آنها نیز مورد سنجش قرار گیرد 


۲- مواد و روش‌ها 

برای شناسایی رفتار دوره‌ای حشکسالی‌های ایران» از آمار بارش ماهانه ۶۱ ایستگاه سینویتیک کشور برای یک 
دوره آماری ۳۱ ساله (۱۹۸۳-۲۰۱۳) استفاده شد. این داده‌ها دارای آماری کامل و قابل اعتماد بودند که اندک خلاء 
آماری آن‌ها با استفاده از همبستگی و مدل‌های رگرسیونی بازسازی شدند. همچنین جهت بررسی صحت و اطمینان 
داده‌های آماری» از آزمون ران تست (1650" 1600) که مورد توصیه مراجع علمی هواشناسی است برای کلیه 
ایستگاه‌ها استفاده شد. نتایج آزمون. همگنی داده‌ها را در سطح بالایی تأیید کرد. توزیع و پراکنش ایستگاه‌های مورد 
مطالعه در شکل (۱) آورده شده است. 
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شکل ۱. نقشه موقعیت و پراکنش ایستگاه‌های مورد مطالعه 


بعد از تشکیل بانک اطلاعاتی مورد نیاز: برای شناسایی فراوانی درجات مختلف خشکسالی‌ها از شاخص 
خحشکسالی موثر" (لال۳) استفاده شد. عانطلذ ۷۷‏ 13 در سال ۱۹۹۹ این شاخص را جهت آشکارسازی 
خشکسالی‌ها و زمان شروع خاتمه و تتش تجمعی آن توسعه داده‌اند. بطور معمول این شاخص, بر خلاف دیگر 
شاحص‌های خحشکسالی؛ با گام زمانی روزانه محاسبه می‌شود (2010 ,۷26 عک 12000167). هر چند که اصول آن 
می‌تواند به داده‌های بارش ماهانه نیز تعمیم داده شود (2008 ,له 6 ۳۵۵۵۵ :2007 ,.1 66 ۷10۲10). طبقه بندی 


شاعص حشکسالی موْثر در جدول (۱) آورده شده استاد 


حدول ۱. طبقه بندی شاخص خشکسالی موثر 


وضعیت نمایه بارش مور 
فوق‌العاده مرطوب بزرگ‌تر یا مساوی از ۲/۵ 
خیلی مرطوب بزرگ‌تر یا مساوی از ۱/۵ 
کمی مرطوب بزرگ‌تر یا مساوی از ۰/۷ 
نرمال کوچک‌تر از ۰/۷ و بزرگتر از ۰/۷- 

خشکسالی ملایم کوچک‌تر یا مساوی ۰/۷- 
خحشکسالی شدید کوچک‌تر با مساوی ۱/۵- 
خشکسالی بسیار شدید کوچک‌تر یا مساوی ۲/۵- 


جهت اطلاع از روش محاسبه شاخص خشکسالی موثر می‌توان به منابع (احتسری و همکاران, ۲۰۰۹: ۱۳0-۱6۵؛ 
کالاماراس و همکاران, ۰۳۰« ۱۲۳-۳) مر اجعه نمود. 

۲-- تبدیل موجک 

برای تحلیل طیفی یک سری زمانی از تبدیلات متنوعی می‌توان سود برد. تبدیل موجک یک تبدیل با آشکارسازی 
قابل تغییر از حوزه زمان به حوزه زمان- بسامد است. با این تبدیل می‌توان اطلاعات بسامدی و زمانی یک سری 
زمانی را بطور همزمان در اختیار داشت. طبق اصل عدم قطعیت هایزنبرگ (اصنهاتهععنا عامتعصتدم) نمی‌توان 
توصیف زمان- بسامد یک سری را به طور دقیق داشت. به بیان ساده بکارگیری تبدیل موجک نشان می‌دهد که یک 
بازه بسامدی غالب در کدام دور زمانی از سری وجود دارد. در محاسبه تبدیل موجک یک سری, از توابع پایه 
استفاده می‌شود که معمولاًتوابع نمایی مختاط هستند. بر این اساس رابطه تبدیل موجک پیوسته را می‌توان به صورت 


ضرب داخلی سری زمانی و یک تابع پایه به فرم زیر نوشت: 


1 ۳۲۵۵۲۷۵ 1۳۲0۱02 
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رابطه (۱) (() ی ا,(ع)) < 4 ند از < (۲,۶) ۱۳۲ (۲:۶) 0۲۷۲۲ 

که در آن: 

رابطه (۲) ات ]سس -() ۷ 
۰ لو 


که در آن ۲ وک به ترتیب پارامترهای انتقال (1:25100) و مقیاس (عط1؟) می‌باشند. انتقال» اطلاعات زمانی 
تبدیل را در بر دارد و پارامتر مقیاس عکس بسامد است. تابع نشان داده شده در رابطه (۲) تابع موجک مادر نامیده 
می‌شود. تمامی نسخه‌های انتقال یافته و مقیاس شده همگی از روی یک تابع اولیه بدست می‌آیند که موجک مادر 
نامیده می‌شود. در تبدیل موجک مقیاس‌های بزرگ متناظر با آشکارسازی بسامدهای پایین و مقیاس‌های کوچک 
متناظر با آشکارسازی بسامدهای بالاست. با توجه به رابطه (۱) که به صورت ضرب داخلی بیان شده است می‌نوان 
این گونه برداشت کرد که تبدیل موجک در حقیقت اندازه‌گیری شباهت بین سری زمانی و توابع پایه (موجک) است. 
به بیان دیگر ضرایب تبدیل موجک پیانگر میزان نزدیکی سری زمانی به موجک در مقیاس مورد نظر است. در این 
پژوهش از نسخه گسسته شده تبدیل موجک با برنامه نویسی در نرم‌افزار متلب استفاده شده است. بعد از انجام تبدیل 
موجک می‌توان بازه‌های بسامدی غالب در سری زمانی را شناسایی کرد. سپس با استفاده از روش‌های متفاوت. 
می‌توان معناداری طیف توان این بازه‌های بسامدی را مورد آزمون قرار داد. در این پژوهش از تبدیل موجک مادر 
مورلت (1۷]01161) استفاده شده است. مورلت مادر حاصلضرب یک موح نمایی مختلط و یک بسته گوسی است که 
به صورت زير نشان داده می‌شود: 

رابطه (۳) 


لته یت 
و هه زد () ,۱۳۷ 


این تابع موجحک پایه است. برای تغییر اندازه موجک بر حسب زمان» از موحک مقیاس استفاده می‌کنيم. موجک 
مقیاس به صورت زیر تعریف شده است: 


ابطه (۶) ۲( - ۶) و۳۳ ۲( -) 
1 ۷ 


که 5 پارامتر مقیاس و 8۲ پارامتر انتقال است. تبدیل موجک (5), 1 یک ضرب داخلی از تابع موجک باسری 

زمانی است. این تبدیل با رابطه زیر داده می‌شود: 
۱ ۱-1 
سس ۷ < () ,۱۷ 
7 

که 2 مقادیر سری زمانی را نشان می‌دهد. علامت ستاره معرف تابع مزدوج مختلط ابا در تحلیل موجک از 

عنیع مرتینز )۱ ۳۸۰ استفاده شده است. همچنین در مطالعات اقلیمی. طیف یس ژزمینه مناسب طیف سرخ است 

( تک 200۳1 رطهعط10 ). این پس زمینه نشان می‌دهد که بیشتر محتوای بسامدی در بسامدهای پایین با 


دوره‌های بازگشت طولانی‌تر قرار دارد. در انتخاب طیف پس زمینه و تعیین سطوح اطمینان از روش تورنسی و کمپو 
و نرم‌افزار مربوطه (۱4۹4۸ 7۱-۸) استفاده گردید. 


۳- تنایج و بحث 
در این پژوهش جهت شناسایی دوره‌های غالب موجود در سری‌های زمانی شاخص خشکسالی موثر از تحلیل 
موجک استفاده شد. بسامد بر حسب دور بر یک ماه یا دوره تناوب بر مبنای ماه سنجیده شده است. با استفاده از 
تحلیل موجک: بازه‌های بسامدی هر دوره زمانی شناسایی گردیدند. به طوریکه توان طیف بازه‌های بسامدی که در 
سطع: ٩۰‏ فرص معاداری فراز اسان و پلترین میران واربانش وا یم مس کردنه استحر ام و مود تنل قراز 
گرفته‌اند. شکل (۲) نتیجه اعمال تبدیل موجک بر سری زمانی شاخص خشکسالی مغثر ایستگاه اردبیل را نشان 


می‌دهد. 


اناع 


دوره (ماه) 


0۴ 10 1 0.1 *10 300 0 زمان ۰ 100 0 
واریانس 
,2 44 1.3 024 00 


شکل ۲. نتیحه اعمال تبدیل موحک بر سری زمانی شاخص خشکسالی مزر (۳01) ایستگاه اردبیل 


می‌دهد که اين سه نوار معنادار هستند. این سه نوار به عنوان سه مولفه اول در جداول هر ناحیه آورده شده است. 
نمودار واریانس نیز نشان می‌دهد که اين سه نوار در سطح ۰ درصد معناداری. واریانس را تبیین می‌کنند. در تحلیل؛ 
این سه نوار تناوبی مورد استفاده قرار گرفته است. همانطور که شکل (۲) نان می‌دهد خطوط پربند دیگری نیز 


۱3۰ جغرافیا و مخاطرات محبطی شماره بیست و پنج 


حسب سال نوشته نشده است. اعداد معرف شماره داده هستند. این فرم از فشردگی بیش از حد شکل در راستای 
افقی جلوگیری کرده و وضوح بیشتری به شکل می‌دهد. جهت ارائه نتایج نهایی پژوهش از ناحبه‌بندی خشکسالی- 
های ایران که توسط نوشیروانی در سال ۱۳۹۳ ارائه شده اننتتت استفاده رن انتون احیهبندی برای همین تعداد 
ایستگاه و برای همین دوره زمانی مورد مطالعه. انجام شده است. نوشیروانی ابتدا شاعص خشکسالی موثر در مقیاس 
روزانه را محاسبه و سپس با استفاده از تحلیل خوشه‌ای اقدام به ناحیه‌بندی حشکسالی‌های ایران نمود. نتایج حاصل 
مکانی خشکسالی‌ها به پنح ناحیه تقسیم‌بندی نمود (شکل ۳. در این پژوهش عناوین زیر برای این نواحی انتخاب 
شدند: ناحیه ۸ ناحیه شمال غرب ناحیه ]: ناحیه کرانه جنوبی دریای خزر ناحیه /): ناحبه نیمه شمالی ایران 


مرکزی, ناحیه 10 ناحیه نیمه جنوبی ایران مرکزی و ناحیه ل: ناحیه جنوب شرقی. 


1 1 ۱ ۱ 


۳ .۰۱ ۸ سس سس ور 
۰ 0 450 300 75150 ۵ 


آَ ۱ ۱ )۱ ۱ آ آ آ آ 
82 80 58 56 54 52 50 48 46 


شکل ۳ ناحیه بندی ایران بر اساس شاخص خشکسالی موثر (نوشیروانی. )۱۳٩۳‏ 


۱-۳- ناحیه شمال غرب 
نتایج تحلیل موجک نشان داد که در این ناحیه فقط در یک ایستگاه بازة بسامدی غالب در کل دوره زمانی یا طول 
سری زمانی وجود دارد. در ایستگاه ارومیه در سطح ۰ درصد معناداری با یک بازه تناوبی ۱۳۲-۰ ماهه با 


مرکزیت تقریبی ۱۳۷۸ ماهه روبرو هستیم که در کل دوره آماری وجود دارد. توان طیف این بازه بالا است. در این 
ایستگاه یک دوره ٩۰-۰‏ ماهه نیز وجود دارد که در بازه زمانی ۱۹۸-۵ خود را نشان می‌دهد. تبدیل موجحک 


سال هفتم 


همچنین نشان داد که ممکن است بازه بسامدی يا بازه تناوبی در طول دوره آماری تغییر کند. در ایستگاه وی در دو 
فورم زمان مبرفار فورمای مشاهاه ی شوه شاد تاو ۱۰ مامت در فاضل رمانی ارات اسان ۱۹۵۳۱۵۱۵۸۵ 
وجود دارد. این بازه تناوبی دارای توان بالایی است. نوسانات در این دوره زمانی به مرکزیت یک دوره ۲۲ ماهه انجام 
گرفته است. این بازه بسامدی در طول دوره آماری و برای سال‌های بعد تغییر کرده و در باز؛ زمانی ۲۰۱۲-۲۰۱۳ به 
تیاه وی اما مت و انس یت و ای ها مان تک ره ۳ مامعن تطف کرت 
در دو ایستگاه تبریز و سنندج دو دورة تناوبی و در ایستگاه‌های سقن زنجان و کرمانشاه یک دورةٌ تناوبی وجود دارد. 
در این ناحیه هیچ ایستگاهی دارای سه دور تناوبی معنادار نیست. در این ناحیه بازُ تناوبی غالبی وجود ندارد و بازه- 
ها کارت زار تایه زا فان تم دهان مفایسه فراوانی توارهای قوش نان ام مها کون فده روز ایم 
ناحیه دوره‌های کوتاه مدت و میان مدت غالب‌اند. اما در دهه سوم این بازه‌های تناوبی تغیبر کرده و به دوره‌های بلند 


مدت نیگن می‌شوند. بنابراین روند تنل گز بازهٌ تناوبی در اين تاحبه مشاهده می‌شود. حدول )۲ دوره‌های 


خوی 


تبریز 


کباتشاه 


حدول ۲ دوره‌های تناویی ناحیه شمال غرب ایران 


بازه اول 
٩2‏ 
۸- ۱۹۸۸ 
۳۹ 
۰- ۱۹۸۵ 
۱-۶ 
۰- ۱۹۸۳ 
۰- ۹۵ 
۱۹۸۳-۲۶ 
۱۰۰-۳۰ 
۵- ۱۹۸۶ 
۹2-۳۰ 
۷- ۱۹۸۶ 
۸- 1۶ 


۱۹۸۳-۲۵ 


بازه دوم 
۰ ۱۳۲ 
۳- ۱۹۸۳ 
۸- ۱۶ 
ار ۱ 
۶-۸ 
۱۹۸-۲ 
۱۰۰-۵۰ 
۱۳ ۱ ۱( 


کاربرد تحلیل موجحک در شناسایی رفتار دوره‌ای خحشکسالی‌های ایران 


دوره تناوب (ماه) 


بازهُ زمانی (سا 


بازهُ زمانی (سا 


بازٌ زمانی (سا 


دورهُ تناوب (ماه) 


بازٌ زمانی (سا 


بازة زمانی (سال) 


دوره تناوب (ماه) 


دورهٌ تناوب (ماه) 


دورهُ تناوب (ماه) 


دوره تناوب (ماه) 


دورهٌ تناوب (ماه) 


ل( 


ل( 


ل‌( 


ل‌( 


ل( 


ل‌( 


۱۲ جغرافیا و مخاطرات محیطی شماره بیست و پنج 


۲-۳- ناحیه کرانه جنوپی دریای خزر 

رفتار ایستگاه‌ها در این ناحیه یکنواخت‌تر است. در ایستگاه اردبیل با یک بازهُ تناوبی که در طول سری زمانی 
پایدار است روبرو هستیم. این باند تناوبی شامل یک دوره ۲۲۰- ۱۲۰ ماهه است که در سطح ٩۰‏ درصد معناداری 
واریانس را تبیین می‌کند. در ایستگاه اردییل دو باند تناویی دیگر هم وجود دارد. اولین بازه یک دورة تناوبی ۵۰-1۰ 
ماهه است که در طول زمانی سال‌های ۱1۹۹۵- ۱۹۸۹ خود را نشان می‌دهد. این بازه تناوبی به آرامی تغییر کرده و در 
یک دوره دو ساله به بازهُ تناویی ۲۶-۲۸ ماهه تبدیل می‌شود. این بازه تناوبی در طی سال‌های ۱۹۹۷-۱4۹۹ بر سری 
زمانی حاکم است. همچنین ایستگاه اردبیل دارای یک دوره تقریباً دو ساله است. در این ایستگاه دوره‌های کوناه 
مدت. میان مدت و بلند مدت مشاهده می‌شود. در اين ناحیه سه ایستگاه اردبیل بابلسر و گرگان دارای سه دور 
تناوبی هستند. بیشترین همخوانی در بین دو ایستگاه بابلسر و گرگان مشاهده می‌شود. دوره غالبی در این ناحیه 
مشاهده نمی‌شود و با قطعیت کمتری می‌توان گفت که روند حاکمیت دوره‌های بلند مدت وجود دارد. دوره‌های میان 
مدت نیز در این ناحیه وجود دارند. در طول سری زمانی وقتی از سال ۱۹۸۳ به ۲۰۱۳ نزدیک می‌شویم تغییرپذیری 


خاصی در دوره‌های تناوبی مشاهده نمی‌شود. جدول (۳) سه بازه تناوبی و دوره زمانی مربوط را نشان می‌دهد. 


جدول ۳. دوره‌های تناوبی ناحیه کرانه‌های جنوبی دریای خزر 


نام ایستگاه بازه اول بازه دوم بازه سوم 
۵۰-۰ ۸- ۲۶ ۱۰-۰ دور تناوب (ماه) 
اردبیل ۱۹۸-۵ ۱۹۹۷-9۹ ۳- ۱۹۸۳ بازة زمانی (سال) 
۱ ۷۰ ۱ سب دورة تناوب (ماه) 
ی ۷- ۱۹۸۳ ۲۰۰۳ ی بازة زمانی (سال) 
۱ ۱۱۰-۳۰ میت ث» دورة تناوب (ماه) 
0 ۱۹۸۷-۲۰ سس سست بازُ زمانی (سال) 
۳ ۵۱-۸۰ سس و دوره تناوب (ماه) 
۸ ۲۰۰ سب تسه بازة زمانی (سال) 
۰- ۶ ۰- ۲ ۲۰-۰ دورة تناوب (ماه) 
۳ ۱۹۹۵-۲۸ ۱۹۹۷-۰۱ ۳ ۲۰۰ بازة زمانی (سال) 
"۳ 1 ۸- 1۶ ۳ دورة تناوب (ماه) 
ن‌ ۸- ۱۹۸۳ ۳ ۱۲ ۲۰۰۹-۷۲ بازُ زمانی (سال) 


۳-۳- ناحیه جنوب شرقی 

نتایج تحلیل موجک نشان داد که در این ناحیه همگنی بیشتری در رفتار تناوبی ایستگاه‌ها وجود دارد. در سه 
ایستگاه فقط با یک دوره تناوبی روبرو هستیم. ایستگاه چابهار دارای دو دور تناوبی است؛ بنابراین تحلیل موجک 
تغییرپذیری دورة تناوبی را فقط در ایستگاه چابهار نشان می‌دهد. این تغییرپذیری از دورة میان مدت به بلند مدت 
است. دو ایستگاه زاهدان و ایرانشهر دارای دوره غالب در کل دوره زمانی مورد مطالعه هستند. برای ایستگاه زاهدان 
این بازه تناوبی یک باند ۱۷۰-۲۱ ماهه است که در سطح ٩۰‏ درصد معناداری می‌باشد. در ایستگاه ایرانشهر این بازة 
تناوبی ۲۲۰- ۱۷۰ ماهه است؛ بنابراین این دو ایستگاه رفتار کاملاً یکسانی را دارند. ایستگاه زابل دارای دوره ۲/۵ ساله 
است. در ایستگاه چابهار دوره ميان مدت و بلند مدت مشاهده می‌شود و در دو ایستگاه دیگر با یک دوره بلند مدت 
روبرو هستیم (جدول 4). 


جدول 4. دوره‌های تناوبی ناحیه جنوب شرقی 


نام ایستگاه بازه اول بازه دوم 

۱۷۰-7 سب دورة تناوب (ماه) 

زاهدان 
۳- ۱۹۸۳ سس بازء زمانی (سال) 
۱۷۰-۰ سسست دورة تناوب (ماه) 
۳- ۱۹۸۳ سس بازء زمانی (سال) 
۲ تحت دورة تناوب (ماه) 

زابل 

۲۰۰۱ سس باز زمانی (سال) 
15_۷7 2-۰ ۱۰۵ دورة تناوب (ماه) 

چابهار 
9- ۱۹۹۵ ۱۹۹۰-۷۲ باز؛ زمانی (سال) 


۳-- ناحیه نيمه حنوبی ایران مرکزی 

در این ناحیه فقط ایستگاه بندرلنگه دارای یک بازهُ تناوبی ۱۳۲-۲۱۰ ماهه معنادار است. این دورهٌ تناوبی در طول 
سری زمانی وجود دارد. بیشترین تعداد بازه تناوبی معنادار متعلق به ایستگاه یزد است. این ایستگاه دارای بازه‌های 
تناوی کزتاه مت استه در این ایستگاه عملا نوارهای تاونی تضیر نمی کنت و بهیشسدگی با بازشتد کین در دوره 
تناوب به چشم نمی‌خورد. مقایسه بازه‌های تناوبی نشان می‌دهد که در ناحیه نیمه‌جنوبی ایران مرکزی بازهُ تناوبی از 
۳- ۲/۵ سال تا بازٌ ۱۳۱۷/۵ سال را شامل می‌شود. در اين ناحیه بک بازهُ تناوبی ٩۶-۸۰‏ ماهه در ایستگاه‌های 
اهواز: آبادان بمء تربت‌حیدریه و مشهد مشاهده می‌شود؛ اما باز؛ زمانی این دوره در ایستگاه‌های یادشده متفاوت 


است. در ایستگاه بوشهر دور تناوبی معنادار مشاهده نشد. رفتار ایستگاه یزد و بندر لنگه با سایر ایستگاه‌های اینن 


۱۹ 


جغرافیا و مخاطرات محیطی 


شماره بیست و پنج 


ناحیه متفاوت است. ایستگاه بندرلنگه بزرگترین بازه تتاوبی رادارد. همچنین در اين ایستگاه یک دوره تقریبا ۳/۵ ساله 


مشاهده می‌شود. ایستگاه یزد دارای یک بازه تناوبی ۶/۵-۵ ساله است که به عنوان مولفه چهارم در بازه زمانی ۲۰۱۳- 
۱ روی داده است. ملفه چهارم فقط در ایستگاه یزد وجود دارد؛ مانند سه ناحیه قبلی» در این ناحیه نیز رفتار 
تناوبی سالیانه مشاهده نمی‌شود. اما به طور کلی می‌توان چنین برداشت کرد که روندی در بازشدگی دوره تناوبی 


وجود دارد و در اواخر طول دوره آماری, نوارهای تناوبی به دورهای بزرگ نزدیک می‌شود. جدول (۵) دوره‌های 


تناوبی و بازهُ زمانی هر ایستگاه را نشان می‌دهد. 


نام ایستگاه 


آبادان 


اهواز 


بوشهر 


جدول ۵. دوره‌های تناوبی ناحیه جنوبی ایران مر کزی 


بازه اول 
- ۵ 
۱۹۸۸-۲۳ 
- ۵ 
۱۹۸۸-۰۱ 
٩-۰‏ 
۱۹۹-۲ 
۹-6۰ 
۱۹۹-۸ 
۲۸-۳۹ 
۷- ۱۹۹۱ 
۳۷/۷ 
۱۹۵۱-0۷ 
۵1-۸۰ 
۱۹۸۲-57 
۱۱۳-۰ 
2-۱ ۱۹۸۵ 
۵5-۸۰ 
۱۹۹۰-۰۶ 
۵1-۸۰ 


۱۹۸-۲ 


بازه دوم 


۱۳۲-۰ 
۱۹۸۳-۳۳ 


۲۰۰۱۲ 


بازه سوم 


بازه چهارم 


دور تناوب (ماه) 
بازة زمانی (سال) 
دور تناوب (ماه) 
بازة زمانی (سال) 


دوره تناوب (ماه) 


بازة زمانی (سال) 


دوره تناوب (ماه) 


بازة زمانی (سال) 


دورهٌ تناوب (ماه) 


بازة زمانی (سال) 


دورهٌ تناوب (ماه) 


بازة زمانی (سال) 


دوره تناوب (ماه) 


بازة زمانی (سال) 
دورهٌ تناوب (ماه) 
بازة زمانی (سال) 


دوره تناوب (ماه) 


بازة زمانی (سال) 


دوره تناوب (ماه) 


بازة زمانی (سال) 


دوره تناوب (ماه) 


بازة زمانی (سال) 


ادامه حدول [ 
نام ایستگاه بازه اول بازه دوم بازه سوم بازه چهارم 
- ۲۶ ۱۸-۸ ۳۲-۰ سب دوره تناوب (ماه) 
۷- ۱۵۸۳ ۹ ۲۰۰۳ ۳ ۲۰ سس بازةُ زمانی (سال) 
۶- ۳/۵ ۲-۵ ۲-۶ سس دوره تناوب (ماه) 
یزد 
۷- ۱۹۸۳ ۸- ۱۵۸۶ ۰ ۱۹۹۰ سك بازةُ زمانی (سال) 


۵-۳- احیه نیمه شمالی ایران مرکزی 

در این ناحیه سه ایستگاه دارای بازةُ بسامدی غالب در کل دوره زمانی مورد مطالعه هستند. همچنین بجز ایستگاه 
کاشان در سایر ایستگاه‌ها حداقل دو دورهٌ تناویی معنادار وجود دارد. ایستگاه‌های شاهرود. شهرکرد و خرم‌آباد دارای 
بازةُ تناوبی معنی‌دار در طول سری هستند. در دو ایستگاه شهرکرد و خرم‌آباد بازهٌ تناوبی کاملاً یکسان است. این بازه 
یک دور تناوبی ۲۲۰- ۲۰۲ ماهه است. در ایستگاه شاهرود این بازم یک دوره ۲۰۸-۲۲۰ ماهه است. بنابراین می- 
توان نتبجه گرفت که این سه ایستگاه دارای بازٌ یکسانی در محدوده ۲۲۰- ۲۱۰ ماهه‌اند. همچنین دومین بازه تناوبی 
ایستگاه شهرکرد و خحرم‌آباد نیز یکسان است. اما دوره زمانی این بازه تناوبی در دو ایستگاه فرق می‌کند. در ایستگاه 
شهرکرد این دوره زمانی سه سال بیشتر طول می‌کشد. در این ناحیه روند بازشدگی یا پهن‌شدگی محسوس نبوده و 
تغییرپدیری پهنای نوار تناوبی در اکثر ایستگاه‌ها مشاهده نمی‌شود. جدول (1) دوره‌های تناوبی در این ناحیه را نشان 
می‌دهد. 


جدول 1. دوره‌های تناوبی ناحیه شمالی ایران مرکزی 


نام ایستگاه بازه اول بازه دوم بازه سوم 

۵- ۶۰ ۸- ۲۶ سب دوره تناوب (ماه) 

سبزوار 
۸- ۱۹۸۳ ۲۰۰۸ سب بازة زمانی (سال) 
۱۹۰-۰ ۵۶-۰ ۲۰۸-۰ دوره تناوب (ماه) 

شاهرود 
۱۵۸۳-۹۸ ۱- ۱۹۸۸ ۳- ۱۹۸۳ بازهُ زمانی (سال) 
۲۰-۸ ۲۰-۰ سس دوره تناوب (ماه) 

سمنان 
۷- ۱۹۹۶ ۹ ۲۰۰ یتست بازٌ زمانی (سال) 
٩1-۰ ۲۰-۸‏ سب دوره تناوب (ماه) 

تهران 
۸- ۱۹۹۶ ۱۵۸۹-۵۹ سس بازةٌ زمانی (سال) 
۶۰-۰ ۳۵ ۰- ۱۰۵ دورة تناوب (ماه) 

فزوین 
2-0۱ ۱۹۸۶ 2-۷ ۱۹۹۶ ۱۹۹۱-2۳ باز زمانی (سال) 
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ادامه حدول ٩‏ 
نام ایستگاه بازه اول بازه دوم بازه سوم 
۷۰-۵ سست 2 دورةٌ تناوب (ماه) 
همدان 
۱۹۸۳ ۳ تیف بازهٌ زمانی (سال) 
۲۰۵-۵ ۰ ۱۱۰ سب دورة تناوب (ماه) 
اراک 
۱- ۱۹۸۳ ۲- ۱۹۹۰ سس بازة زمانی (سال) 
ٍ ۱۲۸-۰ ۰-۰ ۲۰۲ یس دوره تناوب (ماه) 
۱۵۳ ۳- ۱۹۸۳ كت باز؛ زمانی (سال) 
۳۸" ۳1 ء#- 0 سح دوره تناوب (ماه) 
اصفهان 
۸- ۱۹۹۳ ۲۰۰۳-۷۲ سس باز زمانی (سال) 
۱۲۸-۰ ۰ ۲۰۲ سب دوره تناوب (ماه) 
شهرکرد 
- ۱۹۵۸ ۱۹۸۳-۲۳ ِثِ بازهُ زمانی (سال) 
۲۰-۷۶ سس ۳ دوره تناوب (ماه) 
کاشان 
4- ۱۹۸۸ اسب 9 بازة زمانی (سال) 
۶- جمع‌بندی 


نتایج تحلیل موجک نشان داد که خشکسالی‌های ایران دارای رفتار دوره‌ای می‌باشد. رفتار دوره‌ای شاخص 
خشکسالی مثر بسیار هقع بوده و دوره‌های تناوبی از بازه‌های (نوارهای) تناوبی کوتاه مدت تا بازه‌های بلند مدت 
را شامل می‌شود. در اکثر ایستگاه‌ها بازه‌های تناوبی در طول سری زمانی وجود ندارند و فقط در یک بازه زمانی 
معینی رخ می‌دهند. اما در تعدادی از ایستگاه‌ها با بازه‌های تناوبی غالب یعنی بازه‌های تناوبی که در طول سری زمانی 
وجود دارنده روبرو می‌شویم. بنابراین نمی‌توان برای کل ایران بازه تناوبی غالبی را در نظر گرفت. اما در اکشر ایستگاه- 
ها با تغییرپذیری در بازه تناوبی (نوار تناوبی) روبرو هستیم. بجز ناحیه نیمه شمالی ایران مرکزی, در سایر نواحی 
پهنای نوار تناوبی در حال افزایش است. به عبارت دیگر نوار تناوبی به سمت بازشدگی است. در اوایل دوره زمانی با 
نوارهای تناوبی کوتاه مدت روبرو هستیم. در طول سری زمانی و نزدیک شدن به اواخر طول دوره سری» وارهای 
تناوبی به بازه‌های بزرگ نزدیک می‌شوند. بنابراین در کل ایران می‌توان نتیجه گرفت که بازه‌های تناوبی با شیب خیلی 
ملایمی به سمت بازه‌های تناوبی بلندمدت متمایل می‌شود و احتمال رخداد بازگشت خشکسالی‌ها بیشتر می‌شود. ایین 
موضوع می‌تواند بیانگر این باشد که عوامل بزرگ مقیاس نقش موثری در ایجاد دوره‌همای تناوبی بازی می‌کنند و 
نقش عوامل منطقه‌ای رو به کاهش است. برای رسیدن به نتایج دقیق‌تر می‌توان با استفاده از تحلیل طیفی به روش 
فوریه ابتدا دوره‌های تناوبی غالب را شناسایی نمود. سپس با اعمال تبدیل موجک پایداری زمانی اینن دوره‌های 


تناویی را مورد سنجش قرار داد. همچنین برای شناسایی علل این دوره‌ها می‌توان تبدیل فوربه و موجک را بر 
شاحص‌های بزرگ مقیاس مانند الگوهای پیوند از دور اعمال نمود و اثر این شاحص‌ها بر دوره‌های تناوبی شاحص 


خشکسالی موثر را برآورد نمود. 


تشکر و قدردانی 

این پژوهش در قالب طرح پژوهشی با شماره ٩۳۱۲‏ و با استفاده از اعتبارات پژوهشی دانشگاه سیستان و 
بلوچستان انجام شده است. نویسندگان مقاله مراتب سپاس و تشکر ویژه خود را از حوزه معاونت محترم پژوهش و 
فناوری دانشگاه سیستان و بلوچستان که با حمایت مادی و معنوی امکان انجام اين پژوهش را فراهم ساختند اعلام 


می‌دارند. 
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